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Abstract. We consider a general framework for investigating spectral pollu- 
tion in the Galcrkin method. We show how this phenomenon is characterised 
via the existence of particular Weyl sequences which are singular in a suitable 
sense. For a semi-bounded sclfadjoint operator A we identify relative com- 
pactness conditions on a selfadjoint perturbation B ensuring that the limiting 
set of spectral pollution of A and B coincide. Our results show that, under 
perturbation, this limiting set behaves in a similar fashion as the essential 
spectrum. 



Contents 

1. Introduction 1 

2. Limiting spectra 3 

3. Limiting spectra and the behaviour of singular £-Weyl sequences 7 

4. Mapping of the limiting spectra 11 

5. Stability properties of the limiting essential spectrum 14 
References 19 



1. Introduction 

Let A be a self-adjoint operator acting on a separable infinite dimensional Hilbcrt 
space T~L and let A be an isolated eigenvalue of A. For Iclan interval let lz(.A) be 
the spectral projector of A associated to X. The numerical estimation of A whenever 
inf a ess (A) < A < sup<7 ess (^4) and, more generally, when 

Trl ( _ 00 ^)(^) - Trl( A ^ oo) (A) = oo, 

constitutes a serious challenge in applied spectral theory. Indeed, it is well es- 
tablished that classical approaches, such as the Galcrkin method, suffer from varia- 
tional collapse under no further restrictions on the approximating space. This often 
leads to numerical artefacts which do not belong to the spectrum of A, giving rise 
to what is generically called spectral pollution. 

The spectral pollution phenomenon occurs in different practical contexts such as 
Sturm-Liouville operators [1] [32] [31] , perturbations of periodic Schrodinger oper- 
ators [23] and systems underlying elliptic partial differential equations [2] [4] [5] . 
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It is a well-documented difficulty in quantum chemistry and physics, in partic- 
ular regarding relativistic computations [20l [30l El [15]. It also plays a funda- 
mental role in elasticity and magnctohydrodynamics [THl [TUl [HI [3J. In recent 
years this phenomenon has raised a large interest in the mathematical commu- 
nity [231 HH1 HH [HI [131 [HI US]. There are known pollution- free computational 
procedures alternative to the basic Galcrkin method. These include specialised 
variational formulations such as those studied at length in [TBI El El El] as well 
as general methods such as those proposed in [51 (HI EJ ■ 

A natural approach to deal with spectral pollution, is to derive conditions on the 
approximating subspaces guaranteeing a "safe" Galerkin method in a given interval 
of the real line. These conditions were found in [22] on an abstract setting for op- 
erators with particular block-type structures with respect to decompositions of the 
ambient Hilbert space. They turn out to be motivated from techniques in numerical 
analysis [H0 HZ] and computational physics and chemistry (see references in |22j). 

In the present paper we adopt a more general viewpoint than that of |22| . We 
establish an abstract framework for spectral pollution in the Galerkin method and 
then examine its invariance under relatively compact perturbations. Our main 
concern is primarily theoretical and general in nature. Nonetheless, however, we 
include various simple examples which illustrate the many subtleties faced when 
dealing with spectral pollution on a practical setting. 

The technical context of our results can be summarised as follows. Let D(A) be 
the domain of A. Let C = (C n ) n ^ be a sequence of finite dimensional subspaces 
of D(A), dense in the graph norm as n — > oo (Definition [T]Q. Denote by A n the 
compression of A to C n . Denote by a(A,C) the large n limiting set in Hausdorff 
distance of the Galerkin method spectra a(A n ), (Definition [5]). Then <j(A) C 
a(A,C) (Proposition [2]), however in general equality fails to occur in this identity. 
An abstract notion of limiting spectral pollution set can be formulated naturally 
as, 

a pon (A,£) = a(A,£)\a(A). 

As it turns, points in the limiting spectral pollution set behave in a similar fashion as 
points in the essential spectrum (Proposition [3j. Therefore a question arises: what 
sort of conditions on a perturbation B ensure <7 po ii(^4, C) = a po ii(B,C)7 Below we 
establish a theoretical framework in order to address this question. 

Section [2] and [3] are devoted to a characterisation of <r(A,C) in terms of special 
Weyl-type sequences (£-Weyl sequences) and its structural properties. In Defini- 
tion [3J we consider a decomposition of a(A,C) as the union of a limiting essential 
spectrum associated with C, a ass (A, C), and its limiting discrete spectrum coun- 
terpart, adisc(A, £). The former contains both the true essential spectrum a ess (A) 
and a po n(A, C) (Proposition [5]) . 

The purpose of sections [4] and [5] is to find conditions on B ensuring 

(1) <J 0SS (B,C) =a esa (A,C). 



Below we will often consider a slightly more general framework which covers important ap- 
plications such as those involving the finite element method. In this framework we will only 
require that the subspaces C n lie in the domain of the quadratic form associated to A and that 
the sequence C is dense in the form sense. However, in this more general setting we restrict our 
attention to A being semi-bounded. 
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According to our main result (Tlieorcmlll|). when A and B are bounded from below 
and 

(2) (A~a) 1/2 (B -a)- 1/2 -1 

is a compact operator for some a negative enough, ([1]) holds true. Therefore, an 
approximating sequence C will not asymptotically pollute for A in a given interval 
if and only if it does not pollute for B in the same interval. This generalises [22j 
Corollary 2.5]. 

Our present approach consists in adapting to the context of limiting spectra, 
several classical results for the spectrum and essential spectrum. In turns, this 
leads to many unexpected difficulties which we will illustrate on a variety of sim- 
ple examples. In particular, we establish (Theorem [7} a limiting spectra version 
of the spectral mapping theorem allowing to replace the unbounded operator A 
by its (bounded) resolvent (A — a) -1 . Remarkably, this theorem fails in general 
(Remark [4|) for operators which are not semi-bounded. 

2. Limiting spectra 

We will often restrict our attention to A being bounded from below, however we 
do not require this for the moment. Unless otherwise specified, we always assume 
that the subspaces C n are dense in the following precise sense. 

Definition 1 (^4-regular Galerkin sequences). We say that C = (£„), C n C D(A), 
is an A-rcgular Galerkin sequence, or simply an ^4-regular sequence, if for all / € 
T>(A) there exists a sequence of vectors (f n ) with /„ <G C„ such that /„ — > / in the 
graph norm of A, that is: 

(3) ||/n-/ll + IIAf»-4f|| o. 

The orthogonal projection in the scalar product of % onto C n will be denoted by 
7r„ : % — > C n and the compression of A to C n by A n = n n A\c„ '■ C n — > C n . These 
compressions will sometimes be identified with any of their matrix representations. 
On sequences (cc„)„ e N C H of vectors and (£ n ) nS N of subspaces C n C D(A) we 
will often suppress the index and write (x n ) and (C n ) instead, we will denote by 
x n x the fact that x n is weakly convergent to x £ %. When the norm is not 
specified, x„ — > x will denote the fact that \\x n — x\\ — > 0. 

When A is semi-bounded, we may also consider sequences C = {C n ) only in 
the form domain of A. They may approximate the latter but not necessarily the 
operator domain. If A > 0, for instance, this simply means that C n C D(A 1 / 2 ) 
and C is A^-rcgular but not necessarily A- regular. In our notation, for x n G C n , 
Ax„ e D(A 1 / 2 )#, the dual of D(A 1 / 2 ) as subspace of U. Since n n A 1/2 y € T~L* = % 
and , K n A x l 2 y _L g for any y <G H and jeW9 C n , the compression 7r n A \ c n '■ C n — > 
£„ is well defined also in this framework. Moreover, a matrix representation of 
A n can be obtained in the usual manner, via [(A^-^bj ■, A 1 / 2 bu)] d ^i ™ for a given 
orthonormal basis {bj} of C n . We will denote the duality product associated to 
w e D(A^ 2 )* by z i — > (z\w). 

When A is not semi-bounded but its essential spectrum has a gap containing 
a number a, we could as well consider sequences (£„) which are only \A — a] 1 / 2 - 
regular. We have chosen to avoid mentioning quadratic forms for operators which 
are not semi-bounded, because in practical applications (such as those involving 
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the Dirac operator) the domain of \A — a] 1 / 2 does not necessarily coincide with the 
natural domain upon which the quadratic form is defined. 

The limiting spectrum of A relative to the Galerkin sequence £, is the set of all 
limit points, up to subsequences, of the spectra of A n in the large n limit. 

Definition 2 (Limiting spectrum). The limiting spectrum of A relative to C, 
a(A,C), is the set of all A £ K for which there exists Xk £ a(A rik ) such that 
nfc — > oo and Xk — ► A as k — > oo. 

Since all A n are Hermitian endomorphisms, a (A, C) C K. The following lemma 
provides an alternative characterisation of o~(A, C). 

Lemma 1 (£-Weyl sequences). The real number X £ a(A,C) if and only if there 
exists a sequence Xk £ £n k such that \\xk\\ = 1 and 7r„ fc (A — A) Xk — > as k — > oo. 

Proof. According to the definition, A £ a(A,C) if and only if there exists ^ £ I 
and Xk £ C nk with \\xk\\ = 1 such that Xk — > X and 7r nfc (A — Xk) Xk = 0. As 
7r„ fc (A — X) Xk — (Xk — X)xk — > 0, one implication follows immediately. 

On the other hand, let (xk) be as stated. Since the A n arc Hermitian, there 
necessarily exists Xk £ a(A nk ) such that |Afc — A| < \\(A nk — X)xk\ —> 0. Thus 
A £ a(A 7 C) ensuring the complementary implication. □ 

We call (xfe) an £-Weyl sequence for A £ a(A,£), by analogy to the classical 
notion of Weyl sequence [11] . 

Remark 1. Selfadjointness of A n is crucial in Lemma [TJ We illustrate this by 
means of a simple example. Let H = £ 2 (N) and (ej) C H be the canonical or- 
thonormal basis of this space. Let A be the left shift operator defined by the 
condition A : ej i — > e j-i with the convention eo = 0. Let Ck — Span{ei,i < k}. 
For this data an analogous of Lemma Q] is no longer valid. Indeed, if |A| < 1 and 



Xk 



2k A ei ' 



|A| 



then Xk £ C k , \\x k \\ = 1 and 



I _ l\|2 



\Ax k - xx k \ = y ~_ ^ 2k \X\ k ->• 0. 

Therefore any point of the open unit disk is associated with an £-Weyl sequence. 
On the other hand, however, A n is a Jordan block, so a(A n ) = {0} for all n £ N 
and hence necessarily a (A, C) — {0}. o 

The above characterisation of points in the limiting spectrum combined with the 
minimax principle yields the following fundamental statement. 

Proposition 2 (The limiting spectrum and the spectrum). Let C be an A-regular 
Galerkin sequence or, if A > 0, an A 1 / 2 -regular Galerkin sequence. Then, 

(4) a{A)£a{A,C) 
and 

(5) Opoii(A C) := a(A, C) \ a(A) C {f , £+) 
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J — oo for inf a(A) = — oo 
1 inf (Toss (^4) otherwise 

j +oo /or sup ct(A) = +oo 

) supo"ess(^4) otherwise. 

Proof. We start with the general case of an A-rcgular sequence. The classical 
characterisation of the spectrum of sclfadjoint operators ensures that A S o~(A) if 
and only if there is a normalised sequence (y k ) C D(^4) such that || (A — X)yk\\ — > 
(that is (yk) is a Weyl sequence for A). We will now construct an £-Weyl sequence 
from (y k ). According to ([3]), we can find m ) G N 2 such that x k m <E C m , (yk - 

x k n ) — > and (A — \)(yk — i) — > as m — > oo. By virtue of a diagonal process, 
we can extract a subsequence such that ir mk (A — X)x k lk — > 0. Dividing by Ha^JI 
(which does not vanish in the k — > oo limit), gives (U]) as consequence of Lemma [1] 
When A > and C is only A 1//2 -rcgular, a similar proof applies. We explicitly 
take xj^ :— Tr' m yk where ir' m is the orthogonal projection onto C m , for the scalar 
product furnished by the quadratic form associated to A. For all z £ C m with 

Ml = i 

\(z\(A-X)x k m }\ = \(z\(A - X)y k ) + (A + l)(z\y k - x k m )\ 

< \\(A-X)y k l + (A + 1)11^-^11, 

where we have used that (z\(A + l)yk) = (z\ (A + \)-K' m yk) by definition of the 
projection TT' m . Thus 

ll^^-A)^! < \(A-\)y k \\ + (\ + l)\\y k -x k m \\. 

Our assumption that C is J 4 1 / 2 -regular implies that x^ —> yk in D(^4 1 / 2 ) when 
m — > oo. Hence, the desired conclusion is achieved, once again, by a diagonal 
argument. 

The proof of (O is a classical consequence of the minimax principle. It may be 
found, for instance, in J5TJ Theorem 2.1] or [551 Theorem 1.4]. □ 

In [55J Theorem 1.4] the existence of an ^4-regular Galerkin sequence C such 
that a(A 1 C) = [i~,£ + ] is shown. Therefore the inclusion complementary to ^ 
does not hold in general. This is a source of difficulties in applications as there is 
no known systematic procedure able to identify A-regular Galerkin sequences such 
that o~(A) = a(A 7 C). By virtue of ([5]), limiting spectral pollution o~ po \\(A, C) can 
only occur in "gaps" of the essential spectrum. 

Let us now see how a" po n(^4,£) can be characterised in a more precise manner in 
terms of particular £-Weyl sequences. 

Definition 3 (Limiting essential spectrum). We denote by o- css (A, C) the set of all 
A € a(A,C) for which there exists an £-Weyl sequence (xk) as in Lemma [1] with 
the additional property that Xk — 0. 

By analogy to the classical notions, we will call a css (A,C) the limiting essential 
spectrum of A associated to C and the corresponding sequence (xk) a singular 
£-Weyl sequence. 



whe 



f+ — 



Remark 2. From the definition it follows that ct C ss(^4 + K, C) = ct css (^4, C) for any 
sclfadjoint operator K e /C("H). o 
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Definition 4 (Limiting discrete spectrum). The residual set adisc(A, C) = cr(A, C)\ 
cross (A, C), will be called the limiting discrete spectrum of A associated to C. 

We illustrate these definitions by means of various simple examples. 

Example 1 (A a bounded operator). Let H = Span{e^}„ S N where is an or- 
thonormal set of vectors in a given scalar product. Let C n = Spanje^ , . . . , e n _ 1 , /„} 
where /„ = (cos6»)e+ + (sin0)e- for G (0,tt/2). Le10 

A = $>+)<e+|, 

n>l 

that is, A is the orthogonal projector onto Span(e+) and a (A) = <r ess (A) = {0, 1}. 
Thener(A„) = {0, 1, cos 2 9} for all n and a(A, C) = er css (A, C) = {0, 1, cos 2 9}. Here 
x n = e~ is a singular £-Weyl sequence associated to A = 0, x n = e+ is a singular 
£-Weyl sequence associated to A = 1 and x n = /„ is a singular £-Weyl sequence 
associated to A = cos 2 0. o 

Example 2 (A a semi-bounded operator). Let H be as in Example Q] and define 

£„ = Spanje^, . . . , e n _ 1 , e„ }. 

For/±=sin(i) e T± C os(i)e±,let 

A = J2 n2 \fn)(fn\-J2\fn)(fn\ 

which has a 2 x 2 block diagonal representation in the basis (ej). Then cr css (A) = 
{ — 1} and <7disc(^4) = {"- 2 : n € N}. On the other hand 

cr(A„) = i -1, n 2 sin 2 - - cos 2 -, 1, . . . , (n - l) 2 1 , 
y n n J 

where —1 is an eigenvalue of multiplicity n — 1. Therefore 

er css (A, £) = {-1,0} and a disc (A, £) = {n 2 : n E N} . 

The former is a consequence of Proposition G2-(ii) while the latter follows from 
Proposition |5J-(iii) below. 

We can verify directly the validity of the latter as follows. Assume that conversely 
(x k ) was a singular £-Weyl sequence associated with v 2 € <7disc(^4)- Then ir nk (A — 
v 2 )x k — > and Xk — * 0. For m < n k 

p m Tr„ k (A - v 2 )x k = {A- v 2 )p m x k 
where p m = J2i< m \ft)(ft\- Then, on the one hand, 

\\Pn k -iXk - (f+,Xk).f+\\ 2 <\\{A- ^)p nk -ix k f -> 
so that ||p nfc _iXfc|| 2 + || (A — v 2 )p nk ^ix k || 2 — > as k — > oo. On the other hand, 
(A - v 2 )(x k - (e~ k ,x k )e~ k ) = (A - v 2 )p nk _ix k . 

Since | ((A — v 2 )e~, e~) | = | n 2 sin 2 i + cos 2 i — v 2 1 -> |2-^ 2 | > 0, projecting 
each term onto C nk yields (e~ k ,x k ) — > also. But then 1 = \\x k \\ — > 0, which is a 
contradiction, so there are no singular £-Wcyl sequences for v 2 . o 



Here and elsewhere we use the bra-ket notation \f){g\ to denote the linear operator ip i— ► 
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Example 3 (A a strongly indefinite operator) . Let H and C n be as in Example [2] 
Let f,f = ^e+i^e". Let 

^ = EH.c>aii-E^i/« )(/»!■ 

Then a (A) = {±n : n £ N} = CTdi sc (-4)- On the other hand 

a(A, C) = Z, a css {A, C) = {0} and a disc {A, C) = {±n :neN}. 
The proof of the latter is similar to that of the analogous property in Example [2j o 

3. Limiting spectra and the behaviour of singular £-Weyl sequences 

We now examine more closely various basic properties of the limiting spectra 
a {A, C), a ess (A, L) and a d i sc (A, £)■ These properties can be deduced via an analysis 
of the behaviour of different types of £-Weyl sequences. 

Proposition 3 (Limiting essential and discrete spectra and the spectrum). Let C 
be an A-regular Galerkin sequence, or, if A> ; an A 1 / 2 -regular Galerkin sequence. 
Then 

(i) the limiting spectrum <j{A, C) and the limiting essential spectrum a css (A, C) 
are closed subsets of K; 

(ii) moreover a css (A) C a css (A, C) and <7 disc (A,C) C a disc (A). 



Proof. The proof of (i) involves a standard diagonal argument and it is left to the 
reader. For the second statement we need the following auxiliary result which will 
be used repeatedly below. 

Lemma 4. A sequence Xk € Cn k is such that \\xk\\ = 1, Xk x and ir nk (A — 
X)xk — > 0, only when x € Kev(A — A). 

Proof of Lemma^4\ Suppose that (xk) satisfies the hypothesis with C an A-regular 
Galerkin sequence. Let / £ D(^4) and /„ <E C n such that — > f in the norm of 
D(A). Then (ir nk (A — X)xk, fn k ) — > 0. On the other hand, since fk — > / in D(A), 

(ir nk (A - X)x k , f nk ) = (x k , (A - X)f nk ) -s- (x, (A - X)f). 

Thus (x, (A-X)f) = for all / e D(A), so that x G D(A*) = D(A) and (A-X)x = 
as required. 

Suppose now that A > and £ is only A 1 ' 2 -regular. The hypothesis implies 
that (xk) is a bounded sequence in D(A 1 ^ 2 ). Then the proof reduces to the same 
argument, but taking this time f nk in D(A 1//2 ). □ 

We now turn to the proof of (ii) in Proposition [3] The fact that a ess (A) C 
a ess (A, C) is proved similarly to (|4|). It should only be noted that the £-Weyl 
sequence found for A <G er css (^4) additionally satisfies x* — 0. For the inclusion 
Odisc{A,C) C <Tdi S c(^4) note that, if A e a d i SC {A, C), there exists Xk € C Uk such 
that \\xk\\ = 1, Xk — x ^ and 7r„ fc (A — X)xk — > 0. As A ^ a css (A) (by the 
previous part), then either A € o- d isc(A) or A ^ ^(A). By Lemma 21 the latter is 
impossible. □ 

Remark 3. If cr css (v4) = <r css (A, C) then automatically cr d i sc (A) = adisc{A,£) an d 
a(A) = a(A,C). o 
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We will now examine more closely singular £-Weyl sequences associated to points 
A G a css (A,C). 

Proposition 5 (Singular £-Wcyl sequences). Let C be an A-regular Galerkin 
sequence, or, if A > 0, an A 1 / 2 -regular Galerkin sequence. The real number 
A G a css (A, C) if and only if 

(i) either A o~(A) and there exists X k — > A and yk G C„ k such that yf. 
and ir„ k (A - X k )yk = 0; 

(ii) or A £ o~ ess (A) and there exists X k — > A and y k G C Uk such that y k and 
n nk (A - \ k )y k = 0; 

(iii) or A G Odi SC (A) and for any e > 

Rank(l (A _ e . A+£) (A„)) > Rank (l {A} (A)) + 1 
for all n large enough. 

In cases (i) and (iii), A can in some sense be regarded as a point of spectral 
pollution for A relative to C. In case (iii), A G ^{A), but the multiplicity of 
the approximating spectrum o~(A n ) is too large for n large, leading to the wrong 
spectral representation of A in the limit n — > oo. In our definition of the polluted 
spectrum <j po \\(A, C) in ([5]), we have chosen to require that A ^ &(A), following [22] . 
Any A G ct(t4) satisfying (iii) could also be considered as a spurious spectral point. 
However, in case (ii), the singular £-Weyl sequence (yk) behaves like a classical 
singular Weyl sequence. 

Only in cases (i) and (ii) the existence of a singular £-Weyl sequence (y k ) con- 
sisting of exact eigenvectors of A nk such that Tt nk (A — X k )y k = and X k — > A 
is guaranteed. In case (iii) it may occur that all the eigenvectors of A„ k whose 
corresponding eigenvalue converges to A, converge weakly to a non-zero element 
of Kcr(A — A), and that only a linear combination of these eigenvectors converges 
weakly to zero. This can be illustrated by means of a simple example. 

Example 4 (Spectral point satisfying Proposition [5} (iii)). Let H = Span{eo,e„} 
where eo,e^ form an orthonormal basis. Let 

oo oo 
n— 1 n— 1 

Then a(A) = { — 1, 0, 1} and <7 e ss(^4) = { — 1, 1}- The eigenvalue has multiplicity 
one and associated eigenvector eo. Let 

C n = Span{e±, . . . , e±_ ls /±} where /± = 6 ° + Q " 6 " ~ °™ 6 " 



l + (a±) 2 + (aj) 2 



for 



1± 



2(1 T^) 

Then a(A,C) — {0, ±1} = o~ ess (A, C). In this case A n has two eigenvalues ap- 
proaching zero in the large n limit, with corresponding eigenvectors /„ and /~. It 
is readily seen that eo/ y/2 and so only the difference /+ — /~ tends weakly 

to zero. o 
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Proof of Proposition® Let A G o- css (A,C) C cr(A,C). By definition of a(A,C) 
there exists a normalised sequence (yk) such that Tr nk (A — Xk)yk = and Xk — > X. 
The main question is whether one can ensure that yk — i weakly. Up to extraction 
of a subsequence, we may assume that yk — y G Ker(A — A) (by Lemma H}. If 
A ^ cr(A), then Ker(A — A) = {0} and necessarily y = 0, thus (i) follows. 
For the proof of (ii) we require the following auxiliary result. 

Lemma 6. Let V C ~D(A) be a subspace of dimension d > 0, with associated 
orthogonal projector 7Ty . Let e > be such that 

\\irv{A- X)x\\ < e\\x\\ V.t g V. 

There exists N > and a sequence of spaces W n C £„ of dimension d, such that 
for all n> N 

\\ir Wn (A - X)y\\ < 2eVd\\y\\ Vy G W n . 

Proof of Lemma® We firstly assume that C is ^-regular. Let (ej) be a fixed or- 
thonormal basis of V. Then there exists e™ G £„ such that lie™ — e,|L , .. —¥ when 

J II J J II D(-A) 

n — > oo. The Gram matrix G n := ((e™, e"))i<ij<d converges to the d x d identity 
matrix as n — > oo, and therefore, for sufficiently large n, W n := span{e™, j = 
1, d} has dimension d. Now we define an orthonormal basis for W n by 

d 

fe=i 

— 1/2 

Since (G n )kj — > Skj and e™ — > in the graph norm, it is then clear that 

\\ej -f?\\B(A) -> rwoo. 

This shows in particular that ||(7rw„ — 7Ty)A| — > 0, ||7Tyy n — 7ryj — > and hence 
that ||7rw„(^4 - A)/™ - Try (A - A)e J || 0. Let TV > be such that 

||vr w „(A-A)/7|| <2e Vn>N, j = l,...,d. 
For y = Y? ]=1 Vjfj G W„, we get 

d 

\\w Wn (A - X)y\\ < 2eJ2 \w\ < 2eVd\\yl 

3 = 1 

which ensures the desired property. 

When A > and C is only A 1 ' 2 -regular, the proof is the same, using convergence 
in D(A 1 / 2 ) and the fact that ||(7r w „ - tt v )A 1 / 2 \\ -> □ 

The proof of (ii) in Proposition's achieved as follows. Assume that A G o" C ss(^4)- 
For all d G N there exists a subspace Vd C D(A), such that dim Vd = d 2 and 

||(^-%||<^||t/|| Vj/GVrf, 

see for instance [111 Lemma 4.1.4]. According to Lemma [5] and an inductive argu- 
ment, there is a sequence (rid) C N and d 2 -dimensional subspaces Wd C C 7ldl such 
that 

h nd (A-X)y\\ <~||y|| VyeWd. 
This ensures that A nd has at least d 2 eigenvalues in the interval [A — 2/d, A + 2/d] . 
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Let (/™ d )j =1 C Cn d be an orthonormal set of d 2 eigenvectors of A nd , with 

associated eigenvalues {\™ d )f =1 satisfying \\™ d — A| < 2/d. We inductively define 
the following singular £-Weyl sequence for A: 

yi =/r 

V2 = /£ with 1 <S 2 <2 2 such that K^, Vl )\<l/V2 

2/3 =/ 5 ™ 3 3 with 1 < <S 3 < 3 2 such that \(ya,Vi)\ < l/VS for. j = 1,2 

?Ai = fs* with 1 < 5 d < d 2 such that | (y d , Vj ) \ < 1/Vd for j = 1, . . . , d - 1. 
The existence of (5^ is guaranteed by the fact that 

d 2 

VA = l,..,d-l, l = ||y fc || 2 >El(y^/r)l 2 - 

Indeed, there are at most d indices j in the above summation, such that | (yk , f™ d ) \ 2 > 
1/d. Hence, in total, there are at most d(d— 1) indices j such that \(yk, /™ d ) | 2 > 1/d 
for at least one fc = 1, d — 1. Since d(d — 1) < d 2 for d > 1, we deduce that there 
is at least one index j ~: Sd such that \(yk, fj d )\ 2 < 1/d for all k — 1, ...,d — 1. 
By construction \\yd\\ = 1 and \{yi,ijj)\ < l/-^/max(i, j). Thus j/^ — * as fc — > oo, 
ensuring (ii). 

Note that, conversely, if (i) or (ii) holds true, then A G cr css (A,£) by Definition [3] 
Let us now prove that if A G <7 eS s(^4, £) n <7disc(^.), then (hi) holds true. Let 
Xfc G £ nfc be a singular £-Weyl sequence: 7r ni .(A — A)xt — > 0, = 1 and xu 0. 
Let V = Ker(A — A) 7^ {0} and d = dim(V). For n sufficiently large there is a space 
W n C C„ of dimension d such that for all e > 0, there exists N > such that 

Ikn(i4 - A)i/|| < e||v|| Vt/gW„ 

whenever n > N. Let 6>fc = Span{W„ fc , Xk}- Since Xk and W rlfc does not 
increase in dimension in the large k limit, necessarily dim(5fe) = d+1 for all fc large 
enough. For all e > there exists M > such that 

\\7T nk (A-X)y\\<e\\y\\ Vy G S k 

whenever k > M. This ensures that cr(A Ilfc ) n (A — e, A + e) contains at least d+1 
points counting multiplicity and hence the claimed conclusion is achieved. 

It only remains to prove that (hi) implies A G er css (A, C). Each individual eigen- 
vector of A rik might not converge weakly to 0, however there is a linear combination 
of them that does it. We prove this as follows. Let (fj)j=\ be an orthonormal set 
of d + 1 eigenvectors 

A n , /; Kf 1 ; j = i,..., d + i. 

Up to extraction of subsequences we may assume that f!- — ^ fj G Kcr(A — A) for 
all j = 1, . . . ,d + 1. If fj = for some j, then the desired conclusion follows. 
Otherwise, since dimKer(A — A) = d, there exist coefficients (oj) G C d+1 \ {0} such 
that y^j— ? a jfj = 0. Therefore, we may take 
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as singular £-Weyl sequence for A. This completes the proof of Proposition O □ 
4. Mapping of the limiting spectra 

In this section we establish mapping theorems for the different limiting spectra. 
They are a natural generalisation of the analogous well-known result for a (A) and 
a ess (A) (see for example |29[ Section XIII.4]). From now on we assume that A is 
bounded from below and we take C to be an (A — a) ^-regular Galerkin sequence 
with a < inf o~(A). 

Theorem 7 (Mapping of the limiting spectra). Let A be semi-bounded from below 
and let a < inf o~(A). Assume that C is an (A — a) 1 / 2 -regular Galerkin sequence. 
Then 

(6) \eo-(A,£) <^> (X-a)- 1 ea((A~a)-\g) 
and 

(7) Xe<j css (A,£) (X-a)- 1 e<J css ((A-a)-\g) 
where Q = ((A- a) l / 2 C n ) n ^. 

Remark 4. Recall that a sclfadjoint operator A is unbounded (D(A) C Ti) if and 
only if G cr((A — a) -1 ) for one (hence for all) a £ a (A). As it turns out, A is 
unbounded if and only if £ cr csa ((A— a) -1 , g) for one (and hence all) a < mincr(A) 
and [A — a) 1 / 2 -regular sequence C. Formally in Theorem [7] this corresponds to the 
case +oo e cr(A) and (+oo — a) -1 =0. o 

Evidently a result analogous to Theorem [7] can be established when A is semi- 
bounded from above. However, here A is required to be semi-bounded, in order 
to be able to use a square root (A — a) 1 / 2 in the definition of g, and also for a 
more fundamental reason. When a is in a gap of the essential spectrum, it would 
be natural to expect an extension of the above result by considering, for example, 
g = (JA — a| 1 / 2 £„) neN . The following shows that this extension is not possible in 
general. 

Example 5 (Impossibility of extending Theorem [7] for A strongly indefinite) . Let 
% be as in Example [5J Define C n = Span{ef, . . . , e„_ x , cos(0 n ) e+ + sin(0„) e~} 
with e n := tt/4 - A/(2n) for a fixed A G (0, 1). Let 

Then a(A) = {±n : n £ N} = <Tdisc(A). On the other hand 
a(A,C) =o-(A)U{A}, (A,£) = {X} and 

^disc 

Now 

A- 1 =^n- 1 |e+)(e+|-n- 1 |e-)(e-| 

and g = y/\A\C = C. Since A" 1 is compact we have 

a(A-\g)=a(A~ 1 ) and cr eB8 (A-\ Q) = a eaa {A~ l ) = {0}. 

Thus A G cr css (A, C) whereas 1/A ^ a(A^ 1 i g). o 

In fact the following example shows that no general extension of this theorem 
is possible whenever a lies in the convex hull of the essential spectrum, even for 
A G B{H). 
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Example 6 (Impossibility of extending Theorem [7] for a <G Conv{<7 eS s(^4)})- Let 
H = £ 2 (-7r,7r) and Af(x) = sgn(x)/(x) for all / € H. Then a (A) = <r css (A) = 
{±1}. If C is any A-rcgular sequence, then <j(A,£) C [—1, 1]. Fixing a — yields 
(A -a)- 1 = A. Thus also <t(A _1 ) = {±1} and <r((A - a)" 1 , 0) C [-1,1], for any 
A- regular sequence Q. If we construct a sequence £ = (C n ) leading to a spurious 
eigenvalue A € a ess (A,£) n (—1,1), we will always have A -1 ^ a CSB (A~ 1 , Q), no 
matter what Q is. We thus see that Theorem [7] cannot be extended to include a in 
the convex hull of the essential spectrum. o 

Proof of Theorem^ Statement ([5]) will follow immediately from the next result. 

Lemma 8 (Mapping for the spectrum of compressions). Let A be semi-bounded 
from below, let a < inf cr(A) and £„ C D((A — a) 1 / 2 ). Then 

A e (T(w n A\cJ (A — a)" 1 G a{p„(A — a)" 1 \g n ) 

where Q n = (A — a) 1 / 2 £ n and p n is the associated orthogonal projector. 

Proof. Note that A £ a(A n ) if and only if there exists x £ £ n \ {0} such that 

ir n {A - a) 1 / 2 ((A - a)" 1 - (A - a)" 1 ) (A - a) 1 / 2 :! = —^—TT n (A - X)x = 0. 

A — a 

By fixing y = (A — a) 1 / 2 ^ £ (/„ \ {0}, it is readily seen that A £ er(A„) if and only 
if there exist y £ Q n \ {0} such that 

((A - a) 1/2 u, ((A - a)- 1 -{A- a)" 1 ) y} = 

for all u £ £„. Therefore, the statement A € <r(A n ) is equivalent to the existence of 
y £ Gn\{0} such that ((A — a) -1 — (A — a) -1 ) j/ _L Q n which, in turns, is equivalent 
topn^A-o)- 1 -{A-a)- 1 )y = 0. □ 

We now turn to the proof of (JT)). We begin by establishing an alternative char- 
acterisation of the limiting essential spectrum and then we formulate a stability 
result for the limiting spectra with respect to compact perturbations of the regular 
Galerkin sequence. 

Lemma 9 (Alternative characterisation of a ess (A, £)). Let 

T{A):={f{A) : f £ C C (R \ <r css (A), R) } 
T ± {A):={f{A) : /eC c (R\<7 es8 (4),l ± )}. 

Then 

a ess (A,£) = f| a(A + B,£) 

BeF(A) 

(8) = p| <t(A + B,£)= p| a(A + B,£). 

Here C c (f2, R) denotes the set of all real- valued continuous functions of compact 
support in the open set il. Note that J- (A) is a real vector space and J r± (A) are 
cones, all spanned by projectors onto the eigenspaces of A associated with isolated 
eigenvalues of finite multiplicity. At the end of this section it will become clear the 
reason why we highlight the right hand side characterisation in ([5)1. 



WEYL THEOREMS AND SPECTRAL POLLUTION 



13 



Proof of Lemma We only prove the first equality of ((8]) as the proof of the other 
ones follows exactly the same pattern. It is well-known that 

(9) o- esa (A)= p| a(A + B). 

B£F{A) 

Since all the operators in F{A) are of finite rank, then a ess (A + B) = o- CSB (A) for 
all B G J 7 (A). Hence © is equivalent to 

(10) p| a disc {A + B) = 0. 

BeJ r (A) 

From Remark it follows that a css (A + B,C) = a css (A,C) for all B G T(A). 
Therefore a(A + B,C) = o- d isc(A + B,C)U o- css (A, C). Moreover a disc (A + B, L) C 
Cdisc(^4 + B), by Proposition [3l Hence, by (fT0|) . 

P <j disc (A + B,£) C P ^disc (A + B) = 

BeT(A) BeJ^iA) 

and the result is proved. □ 

Lemma 10. Let T = T* be such that \\T\\ < oo and let C be a T-regular sequence. 
Let A e K(U). If -1 g a(K), then 

a css (T,C) = <t css (T,(1 + K)C) and a disc (T, C) = o- disc (T, (1 + K)C). 

Proof of Lemma \1(A We firstly prove that 

(11) a css {T,C)\ (T) 

^"css 

(T,(l + K)£)\ 

^"disc 

Since 

(12) C = (1 + TiT)— x (l + K)C = (1 - A(l + J0 _1 )(1 + K)C, 

it suffices to show that the left hand side of (fTTj) is contained in the right hand side. 
Let A G cr css (T, C) \ o"disc(7 1 ). If A £ a css (T), a direct application of Proposition [3} 
(ii) ensures that A lies also in the right hand side of (jlip . so we can assume that 
A g" <t(A). According to Proposition [5} (i), there exists A& — > A and x k G £ nfc such 
that ||xfc|| = 1, Xk — 1 and 7r„ fc (T — Afc)xfc = 0. For all v k G £n fe , and hence for all 
w k = (1 + A> fc G (1 + K)C nh , we have 

0= ((T- \k)xk,v k ) = ((l + K*)- 1 ^- \ k )x k ,(l+K)vk) 

= ((l + A'^-^T-AfeK,^). 

Let q k be the orthogonal projection onto (1 + K)C nk . Then 

q k (l+K*)- 1 (T-\ k )x k =0. 

Now (1 + K*)- 1 = 1 - K where A = A*(l + A*)" 1 G K.(H). Hence 

q k (l-K){T-X)x k ^0. 

But, since ||Tj| < oo and x k — 0, A(T — A)xfc — > 0, so that also gfc(T — X)x k — > 0. 
Thus q k (T — X)y k — > for j/^ = (1 + K)x k 0. By re-normalising y k in the obvious 
manner, we obtain a singular £-Weyl sequence for A G cr(T, (1 + A')£), ensuring 

To complete the proof of the first identity in the conclusion of the lemma, suppose 
that A G <J CSS (T,C) n a disc (T). For any fx £ A let f = T + (jj, - A)l (A _ e>A+e) (T) 
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where e > is sufficiently small. Then A G a css (T, C) \ o-<n$c(T). By virtue of (fTTj) 
and Remark El A € er 0SS (T, (1 + A')£) = cr ess (T, (1 + A")£) as needed. 

We now show the second identity in the conclusion of the lemma. By virtue of 
(|12p and the first identity which we just proved, it is enough to verify 

<7 disc (T,£) C a(T,(l + K)C). 

This, in turns, follows from Proposition [3]- (ii) and (HJ, since 

(T,£)c (T) and er(T) C cr(T, (1 + A)£) 

taking into account that (1 + A)£ is a T-regular sequence. □ 

We now complete the proof of Theorem [7] by showing (0. Let A G <r css (A,C). 
By virtue of Lemma [9] this is equivalent to the statement 

\/BeF+(A), \e<t(A + B,£). 

Since B > and a < min[cr(A + B)], according to (j6|) the latter is equivalent to 

VBeJ+(A), (A -a)" 1 e a((A + B-a)-\g B ) 

where £ B = {A + B - a)- 1 / 2 ^. Since -B has finite rank and is therefore compact, 
Lemma [10] ensures that the above in turns is equivalent to 

VB e F + {A), (A - a)" 1 G a((A + B- a)" 1 , So). 

Note that ^ + B — a) 1 / 2 {A — a) -1 / 2 ) as the corresponding operator is 

an invertible function of A. Now (A + B — a) -1 = (A — a) -1 + B, where B = 
-(A~a)- 1 B(A+B-a)~ 1 runs over all of T~ {(A - a) -1 ) as B runs over all T + (A) 
and conversely. For the latter note that / G J 7+ (A) if and only if— /((• — a) _1 )G 
J-~ {{A — a) -1 )- Thus, once again by Lemma 03 A G cr 0SS (A,£) is equivalent to 

(X-a)- 1 ea css ((A-a)-\G). 

This completes the proof of Theorem [7] □ 

Remark 5. The above proof mimics the proof of the classical Mapping Theorem 
for the essential spectrum which can be deduced from the characterisation 

<j css (A)= p| a (A + B), 

BeK,{U) 

see, e.g., [15]. o 

5. Stability properties of the limiting essential spectrum 

In this final section we present the main contribution of this paper. It strongly 
depends on the validity of Theorem [7] 

Theorem 11 (Weyl-type stability theorem for the limiting spectra). Let A and 
B be two selfadjoint operators which are bounded below. Assume that for some 
a < m£{a(A),<T(B)}, 

(13) D((B-a) 1 / 2 ) = D((yl-a) 1 / 2 ) 
and 

(14) (A - a) 1/2 ((B - a y 1/2 - (A - a)- 1 ' 2 ) G JC(H). 
Then 

a css (A,C) = cr css (B,C) 
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for all sequences C = (£„) which are simultaneously (A — a) 1 / 2 -regular and (B — 
a) 1 / 2 -regular. 



Under Assumption (|13[) . (|14[) is equivalent to the same condition with the roles 
of A and B reversed: 

(15) (B - a) 1/2 ((A - a)- 1 ' 2 - (B - a)- 1 ' 2 ) G JC(U). 



Note however that ([13]) and ([14]) do not imply necessarily that an A-regular sequence 
is also B-regular. For this it is enough to consider an example where D(A) 7^ D(B). 
Let A = <9 4 with domain 

D{A) = B 4 (0, 1) n {u(0) = u{l) = 0, u"(0) = u"(l) = 0} C L 2 (0, 1) = H. 

Let B = + |1)(1|(1 - dl) with domain 

D(B) = B 4 (0, 1) n |u(0) = u(l) = 0, u"(0) = u"(l) = J u(x)dx| 7^ D(A). 

Then A 1 / 2 = -<9 2 and B 1 / 2 = -d 2 x + |1)(1| both with domain 

D(A 1/2 ) = D(B 1/2 ) = B 2 (0, 1) n {u(0) = u(l) = 0}. 

Since A 1/2 (B~ 1//2 - A~ 1/2 ) = |1)(1|B~ 1/2 is a rank-one operator, A and B satisfy 
the hypotheses of Theorem [TT1 but clearly A-regular sequences C n C D(A) \ D(B) 
are not B-regular. 

Remark 6. The KLMN theorem [2§] ensures that if B — A is a densely defined 
symmetric A- form-bounded operator with bound less than 1, then (|13[) holds for a 
sufficiently negative. 

The following example from [22] shows that Theorem [11] cannot be easily gener- 
alised to operators which are not semi-bounded. 

Example 7 (Relatively compact perturbations of the Dirac operator) . Let A = D° 
and B = D° + V where D° denotes the free Dirac operator with unit mass [55] 
and V £ C^°(R 3 ) is a smooth non- negative function of compact support. The 
ambient Hilbcrt space here is H — i 2 (M 3 ,C 4 ). Under the additional assumption 
that supU = lV| ia o(R3) < 1, h is guaranteed that ^ &(B). Furthermore it can 
be verified that 

D (V| 1/2 ) = D (|B| 1/2 ) = # 1/2 (R 3 , C 4 ) 

and that 

| A |i/ 2 (|5|-i/2 _ l^l" 1 / 2 ) £JC(H). 

As a consequence of (22] Theorem 2.7], it is known that there exists a B-regular 
Galcrkin sequence C = (£„) such that 

(16) a css (B,C) D [0 ; supV]. 

These Galerkin spaces comprise upper and lower spinors, meaning that 



for suitable (/"), (<?") C L 2 (M. 3 ,C 2 ). This basis is known to be free of pollution if 
the external field V = 0, that is 

a(A, C) = a(D°) = (-00, -1] U [1, 00) = a css (A, C). 
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Hence a ess (A,C) ^ <j e8S (B,£) so Theorem [TT1 fails for operators which are strongly 
indefinite. o 

Proof of Theorem [771 Denote by K the operator on the left side of (fT4f . Then 

(17) (B - a)- 1 -{A- a)- 1 = {A - a)- l ' 2 K{B - a)- 1/2 + (A- a)- x K e fC{H). 
Let Q := (A - a) l / 2 C. According to ©, 

Xea ess (A,C) (X-a)- 1 eo- css ((A-a)-\g) . 

By Remark [21 

a css {{B - a)-\g) = a css ({A - a)" 1 , S) . 

Let Q 1 = {B- a) l / 2 C. Then 

g = (A - a)^ 2 C = (A - a)^ 2 (B - a)- 1 ' 2 g l = (1 + K)g> '. 

Note that K = {A — a) 1 ^ 2 (B — a) -1 / 2 — 1 and —1 ^ cr(K) as a consequence of the 
fact that £ a((A - a)^ 2 (B - a)~ 1/2 ) by (JT3J) - According to Lemma[10l 

(18) a css {{B - a)-\g) = a css ((B - o)" 1 ,^) 

The conclusion follows by applying Theorem [7] again, this time to operator B. □ 

Corollary 12. Let A and B be two bounded-below selfadjoint operators such that 
(|13p holds true for some a < mi{a(A),a(B)}. Assume thatC := B — A is a densely 
defined symmetric operator such that 

(19) C eB(D((B-af),H) 
and 

(20) {A-a)- a C{B-a)- p e K(H) 
for some < a, f3 < 1 with a + [3 < 1. Then 

cr oss (^4, C) = a css (B,C) 

for all sequences L = (C n ) which are simultaneously (A — a) 1 / 2 -regular and (B — 
a) 1 / 2 -regular. 

Remark 7. Let A be a given bounded-below selfadjoint operator and assume that 
A has a gap (a, b) in its essential spectrum in the following precise sense, 

<j css (A)D(a,b) = 0, tr (l^^A)) = tr (l (6j0o) (^)) = +oo. 

Let LI := 1( Cj00 )(t4) where a < c < &. Results shown in [55] ensure that, when the 
Galcrkin spaces C n are compatible with the decomposition W = ILH©(1 — 11)% (i.e. 
when II and 7r„ commute for all n), there is no pollution in the gap: a css (A, C) fl 
(a, 6) = 0. According to [551 Corollary 2.5], when 

(21) {B-a^CiA-a)- 1 ' 2 € 1C(H), 

then a ess (B,£) = as well. 

In this respect, Theorem [TT] can be seen as a generalisation of these results. 
Although condition (|20|) is stronger than pip , the statement guarantees that the 
whole polluted spectrum will not move irrespectively of the (A — a) 1 / 2 -regular 
Galerkin family C and not only for those satisfying [II, 7r n ] = for all n. o 
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Example 8 (Periodic Schrodinger operators). Let A = —A + V pci - where V pcl is a 
periodic potential with respect to some fixed lattice KcR' i (for instance 7Z = Z 3 ). 
Let C = W(x) be a perturbation. Assume that 

[ p = 2 if d < 3 
V pei e Lf oc (R d ) where < p>2 if d = 4 
[ p = d/2 if d > 5 

and that 

W G L 9 (M d ) n Lf oc (K d ) + L™(R d ) 
for max(d/2, 1) < q < oo. Then (|20[) holds true for suitable a, ft and a, and 
therefore 

(22) CTcss (-A + y pcr + W, C) = <j css (-A + T/ por , £) 

for all A-regular Galerkin sequence C. See |22[ Section 2.3.1]. 

A Galerkin sequence C which does lead to any pollution in a given gap, can 
be found by localised Wannier functions, [221E]. In practice, these functions can 
only be calculated numerically, so it is natural to ask what would be the polluted 
spectrum when they are known only approximately. According to (|22[) . the polluted 
spectrum will not increase in size more than that of the unperturbed operator 

-A + Fper- O 

Example 9 (Optimality of the constants in Corollary [T2")) ■ Let H, £, and 
be as in Example [2j Let 

A = J2n l \f+){ti\ + J2\fn){fn\ and B = ^n r |e+>(e+|+^|e-)(e-|. 

n n 

The matrix representation of A and B in the basis is made out of 2 x 2 blocks 
placed along the diagonal. More precisely A = diagL4„], B = diag[£?„] and C = 
diag[C„]; where 

A - P ( nt °\ R R — ( nT ° 

o i / — \ 1 



and C n — A n — B n for 



cos — sin i 



Fix a = and let L — A a CB 13 . The matrix representation of L in the basis et 
is L = diag[L„] where we can calculate explicitly the entries as 

(L n ) n = -n-^-at+r CQS 2 I + n -?r ^2 I 

n n 
- n- r ^ sin 2 - + n-^- aW cos 2 - 



n 



n < a ~^ - n- ae 



, s 1 1 ( 

(in) 12 = cos - sin — 

n n \ 

(L n ) 21 = cos - sin - ( n - fir - at+t - n -^- ai+r ~ n^ r + n~ r ^ 
n n \ 



n 

{L n ) 2 2 = sin 7 1 
n 



2± r„-Ao-i) 
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Therefore L is compact, given the following 

1 = 2, 0</3,a<l, 0<r<2, 

(23) 1 
-/3r-2a + 2<0, a > 1/2, /3 > 1 - -. 

r 

On the other hand, for I = 2, 

aes S (A,£) ={1,2} and a 0SS (B, £) = {1}. 

This example suggests that condition (|2T)|) in Corollary is quasi-optimal for the 
stated range of (3 and a as illustrated by Figure [1] Note however that in this 
example (fT3")) is only satisfied when r = I. o 




Figure 1. The region in green ocolour for the parameters /3 and 
a is covered by the conditions of Corollary [T^l If A and i? satisfy 
(|20[) for (/3, a) in this region, then the limiting essential spectrum 
is preserved. The region in red shows the parameters f3 and a in 
condition (|23| of Example [9] The region in blue is generated by 
exchanging the roles of (3 and a. It is not enough for A and B to 
satisfy (|20|) for (/3, a) in these two regions, to guarantee preserva- 
tion of the limiting essential spectrum. 



Proof of Corollary 1 121 Assume firstly that < a < 1/2. The proof reduces to 
showing that the operator K defined by expression (|T4)) is compact. Let L = 
(A - a)~ a C{B - a)- fi be the operator given by ([2"0]l. Since (3 < 1, we have D(B - 
a) C D(S - a)^ 3 , [Til Theorem 4.3.4]. Then, by C D((A - a) Q ) and 

Ca; = (A - a) a L(B - a) /3 a; for all x G D(B - a). By virtue of 0, 

(i4 - a) 1/2 (A - a + 8 )- x C(B - a + s)" 1 6 iC(W) 

for all s > 0. Moreover 

(A-a) 1/2 ((A-a + s)- 1 -(B-a + s)- 1 )a- = (^-a) 1/2 ( J 4-a + s)- 1 C(S-a + s)- 1 a; 
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for all x € H, as this identity is satisfied in a dense subspace of H. Thus 

1 /*°° A 

K = / (A — a) 1/2 M - a + s)~ 1 C{B -a + s)" 1 -^ 

= -- !°° \{A- ay/ 2 (A -a + S )-\A - a) a ) L {{B - af {B - a + s)- 1 } 
^ Jo 1 J V s 

Both terms in brackets multiplying L arc bounded operators, then the integrand in 
the second expression is also a compact operator. Moreover, the integral converges 
in the Bochner sense as its norm is 0(s l3+a ~ 2 ) for s^oo and 0{s~ 1 / 2 ) for s-)0. 
Thus K e fCi^hL) in this case and Theorem ITT1 implies the desired conclusion. 
Now suppose that 1/2 < a < 1, so that < /3 < 1/2. Since 

B(A - a) a c D(A - a) 1/2 = B(B - a) 1 ' 2 c B(B - af , 

then C € #(D(,4 - a) a ,H). Hence the operator (B - a)-^C{A - a)~ Q is bounded 
and (B — a)^^C(A — a)~ a x = L*x for all x € H. The proof is then completed by 
exchanging the roles of A and B. □ 
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